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The presence of (approximate) conservation laws can prohibit the fast relaxation of interacting 
many-particle quantum systems. We investigate this physics by studying the center-of-mass oscil¬ 
lations of two species of fermionic ultracold atoms in a harmonic trap. If their trap frequencies 
are equal, a dynamical symmetry (spectrum generating algebra), closely related to Kohn’s theorem, 
prohibits the relaxation of center-of-mass oscillations. A small detuning 5uj of the trap frequencies 
for the two species breaks the dynamical symmetry and ultimately leads to a damping of dipole 
oscillations driven by inter-species interactions. Using memory-matrix methods, we calculate the 
relaxation as a function of frequency difference, particle number, temperature and strength of inter¬ 
species interactions. When interactions dominate, there is almost perfect drag between the two 
species and the dynamical symmetry is approximately restored. The drag can either arise from 
Hartree potentials or from friction. In the latter case (hydrodynamic limit), the center-of-mass 
oscillations decay with a tiny rate, l/r oc {6ujY/V, where T is a single particle scattering rate. 

PACS numbers: 67.85.-d, 67.85.Lm, 61.20.Lc, 67.10.Jn 


How does an interacting many-body quantum system 
reach thermal equilibrium? While often a few scattering 
processes are sufficient to establish locally an approxi¬ 
mate equilibrium state, in some cases the presence of 
conservation laws prohibits equilibration. In one dimen¬ 
sion (ID), for example, integrable quantum systems like 
the spin 1/2 Heisenberg model or the fermionic Hub¬ 
bard model possess an infinite number of conservation 
laws. Due to their presence, the system cannot relax to 
a simple thermal state described by just a few parame¬ 
ters like temperature or chemical potential. Instead, only 
an equilibration to a generalized Gibbs ensemble (GGE) 
[1] is expected where for each conservation law a new 
Lagrange parameter is needed to describe the long-time 
steady state. 

Real experimental systems are, however, often only ap¬ 
proximately described by integrable models. As a con¬ 
sequence the corresponding conservation laws are only 
approximately valid. For classical systems with a finite 
number of degrees of freedom, the famous KAM theorem 
[2H4] states that even in such a situation many proper¬ 
ties of the integrable point can survive. The situation for 
interacting many-particle quantum systems is less clear. 
Generally it is, however, expected that due to integrabil- 
ity breaking terms the system can relax to a thermal state 
but the relaxation is slow and governed by the slow re¬ 
laxation of the approximate conservation laws. A similar 
question arises in transport studies: Integrable systems 
like the ID Heisenberg model are characterized by infinite 
(heat-) conductivities even at finite temperature [SHZ]- 
In real materials, however, small integrability breaking 
terms can render the conductivity finite. This has mo¬ 
tivated early studies of the role of integrability breaking 
terms for transport properties 019]. 

Ultracold atoms provide new opportunities to inves¬ 
tigate the question of equilibration and the role of (ap¬ 
proximate) symmetries. For example, in a famous ex¬ 


periment termed “Quantum Newton’s Cradle” uni it was 
shown that the breathing mode of a ID Bose liquid in a 
harmonic trap does not relax on experimentally relevant 
time scales. While in this case the harmonic traps nom¬ 
inally break integrability this apparently has little effect 
on the experiments. 

In this paper, we study equilibration in the presence 
of an approximate symmetry in a model which is (i) ide¬ 
ally suited for experimental studies and (ii) conceptually 
simple due to the presence of only a single symmetry - 
instead of infinitely many. We study the center-of-mass 
(COM) oscillations of atoms in a harmonic trapping po¬ 
tential. If all atoms have the same mass and same trap¬ 
ping potential, then the COM oscillation never decays 
and its frequency is exactly given by the non-interacting 
result m- A closely related results is Kohn’s theorem 
[12] stating that cyclotron resonances of electrons in a 
Galileian invariant system are not affected by interac¬ 
tions. Mathematically, this can be traced back to fact 
that the total momentum P, the center-of-mass R and 
the interacting many-particle Hamiltionan in the pres¬ 
ence of a trapping potential ^Vbr^ form a closed algebra 
(a so-called speetrum generating algebra) given by 

[R\P^]=ih5ii, \P\H] = -ihVoNR\ [R\H]=if^P^ 

Nm 

( 1 ) 

This algebra implies that the COM motion completely 
separates from all many-particle excitations in the trap 
even in the presence of a time-dependent trapping po¬ 
tential ^Vo{t){r — 'ro(t))^. Furthermore, the nonlocal (!) 
operator Q = P^/(2m) + is a conservation law, 

[Q^H] = 0. We will study how these symmetries break 
down when two species of atoms with slightly different 
masses or slightly different trapping potentials are con¬ 
sidered. Such a case has recently been studied by the 
Salomon group [13] using mixtures of ^Li and ^Li. This 
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work investigated, however, mainly the role of superflu¬ 
idity on the COM oscillations in this system. 

The case of a fermionic mixture has been studied the¬ 
oretically by Chiacchiera, Macri and Trombettoni m- 
While the methods used by the authors are similar to 
the one used in our study (projection on the dynamics of 
slow modes), their paper mainly focuses on counting the 
number of relevant modes and contains little information 
on the question discussed in this paper, especially on the 
behavior of the damping rate as function of the trap- 
frequency difference, population difference, inter-species 
scattering rate, and temperature. Furthermore, it only 
considers the classical high-temperature limit where ef¬ 
fects of Pauli blocking can be ignored. 

An alternative option to perturb the dynamical sym¬ 
metry of COM oscillations is to consider corrections to 
the confining harmonic potential, e.g., by adding a 
term. Such a situation has very recently be investigated 
for a quasi one-dimensional setup in the hydrodynamic 
limit by Iqbal, Levchenko and Khodas [15] using the 
Navier-Stokes equation. Similar to the case discussed 
in this paper, they obtain a long-lived mode where the 
decay rates are controlled by the strength of the anhar- 
monic terms. 

A spectrum generating algebra also characterizes ap¬ 
proximately the breathing mode (monopole oscillations) 
of atoms in a harmonic trap in two dimensions [16] and of 
a unitary gas in arbitrary dimensions m- In this case, 
the shift of the resonance frequency due to deviations 
from the unitary limit has been calculated in one dimen¬ 
sion in Ref. naiin] while the two-dimensional case was 
studied in [2Qll22] and investigated experimentally in m- 
Within our study, we will be mainly interested to study 
the relaxation rate rather than the frequency shift. 

The relative motion of two species of atoms is mainly 
controlled by their mutual interactions. This problem, 
often described by the term spin drag^ has been investi¬ 
gated both in the context of electrons in solids [24| [25] 
and also for ultracold atoms, see e.g. [261127]. 

In the following, we will first introduce the model and 
our analytical approach, identify three important physi¬ 
cal regimes (ballistic, frictionless drag, and friction domi¬ 
nated drag), and, finally, quantitatively predict how these 
regimes determine properties both for the real-time evo¬ 
lution and for the response as function of the frequency. 


I. MODEL AND METHOD 


In this article we study two species of ultracold 
fermions with creation operators T|(r) and ^^(r) cap¬ 
tured each in a perfectly harmonic trap in three dimen¬ 


sions. The system is described by 




H 


( 12 ). 
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( 2 ) 


^int ~ 


2777/red 


J d^r T|(r)T2(r)T2(r)Ti(r) 


Here, a shift of the position of the potential minimum 
r^(t) can be used to excite dipolar oscillations. In gen¬ 
eral, the two fermion species may have different masses 
rui and feel different trap potentials with respective trap 
frequencies uji = 00 + ^^ 002 =00 — ^. describes the 

intra-species interaction which we do not specify here as 
it does not influence our results in any qualitative way. 
Furthermore, for spinless fermions safely be ne¬ 

glected. As we will show, all relaxation arises from the 
inter-species interaction which we parametrize by the s- 
wave scattering length a with rUred = 
being the reduced mass (note that we use a pseudopo¬ 
tential to describe the scattering, see, e.g.. Ref. [28]). For 
= 0 the COM oscillations do not decay (see below). 
We are therefore mainly interested in the limit Suj u, 
where a slow decay of the oscillations can be expected. 
Experimentally, this can, for example, be realized by us¬ 
ing two isotopes with slightly different mass, m ± 
but identical trapping potential. In this case ^ 
Alternatively, one can use two hyperfine states of the 
same atom in combination with a spin-dependent poten¬ 
tial [28] . The latter setup has the advantage that one can 
directly tune the parameter 

Our theoretical approach is based on the idea that for 
^ ^ 1 the dynamics is governed by an approximate dy¬ 
namical symmetry which prohibits a fast relaxation of 
the COM oscillations. Furthermore, in the limit of van¬ 
ishing inter-species interactions, a ^ 0, also the COM 
motion of each atomic species separately decouples. Our 
central goal is to derive an effective, hydrodynamic de¬ 
scription of the slowly relaxing modes. We will therefore 
focus on the dynamics in the operator space spanned by 
the center-of-mass coordinates and the total momen¬ 
tum Pi of each of the two species defined by 


R. = Ty‘dVv|/t(r)rv|/,(r); (3) 

Pi = [d^r «'J(r) {-ihV) 


where Ni is the number of particles of type i = 1,2. 

For weak excitations of the system, it is sufficient to 
study linear response within the Kubo formalism. The 
main goal is thereby to calculate the matrix of retarded 
susceptibilities 

• poo 

Xmn(w) = ^ J e*“‘([A„(t),^„(0)])eq. (4) 

where (•)eq. denotes the expectation value in equilibrium 
for v°{t) = 0 Vt and An = ,P^). As for a 
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spherical potential the x, y and 2 ; components do not 
mix within linear response, we can focus on the x co¬ 
ordinate only. Xmn allows to calculate all experiments 
where the COM oscillations are excited by a shift r^(t) 
of the potential and where the COM and/or the average 
momenta of the particles are observed. 

To calculate Xmn(^) we use the so-called memory ma¬ 
trix formalism |29l43T] . The memory matrix is a matrix 
of relaxation rates of slow variables, which we evaluate 
perturbatively in the strength of the inter-species inter¬ 
action. This formalism has the advantages that (i) it is 
easy to evaluate - without the need to solve the type 
of integral equations needed for Boltzmann approaches 
or when vertex corrections are taken into account within 
the Kubo formalism, (ii) it nevertheless automatically in¬ 
cludes the effect of vertex corrections, which are essential 
to describe momentum conservation, which is also gov¬ 
erning the COM oscillations [32], (iii) it is accurate in 
cases where there is a separation of time scales and all 
slow modes are included in the memory matrix, (iv) it 
can be used to treat complicated situations like the ex¬ 
pansion around a fully interacting integrable system 
and has recently been used to calculate transport prop¬ 
erties of exotic non-Fermi liquids j33h35] (v) in the case 
considered here, where we effectively expand around the 
non-interacting limit, it is equivalent to a solution of the 
Boltzmann equation by projection onto the slow modes 
misni. In Ref. |37| we have argued that the formalism 
gives always a lower limit for conductivities. The situa¬ 
tion investigated here is, however, more complicated com¬ 
pared to the case considered in Ref.|37jas we are studying 
here effects at finite frequency in a system which is not 
translationally invariant. This leads to extra dephasing 
effects discussed in detail in Appendices |B 4| and [C] 

We refer to Appendix for a brief review of the 
memory-matrix method. It allows to express the matrix 
Xmn{^) of retarded susceptibilities (cf. Eqs.( A4|A^ ), 


x{^) = — cj (cj — Co 


( 5 ) 


in terms of an equal-time correlation matrix Co, a con¬ 
stant matrix Q and a frequency-dependent matrix-valued 
complex function S(co’). The latter two matrices have 
a similar role as the self-energy: they describe directly 
the shift of frequencies and the damping of oscillations. 
They have the advantage that they can be evaluated di¬ 
rectly in perturbation theory, without the need to resum 
an infinite series of diagrams. More precisely, the latter 
statement holds in the case when all slow modes have 
been included in the set of observables A^. We will use 
An = (Rf, Rf 5 P 2 ) slow modes, which is suffi¬ 

cient to describe the regime where interactions dominate. 


As we discuss in detail in section [B^ of the appendix, in 
the limit of vanishing interactions an infinite set of fur¬ 
ther slow modes exists, which have to be included to 
describe details of the dephasing of oscillations for very 
weak interactions (ballistic regime) studied in detail in 
Appendix [O but not captured for the above choice of 


In the following our goal will be to calculate for weak 
interactions the frequencies and decay rates of the center- 
of-mass oscillations. In appendix we evaluate the ma¬ 
trices E(u;), and Co in local density approximatio n for 
weak int eractions. To linear order in a, using Eqs. (B 6 ) 
and (B9), we find Vt = with 





(6) 


and 


z= i^M20J2 



( 7 ) 


Here, Mi = Niirii is the total mass of the fermions of 
species i and 7 has the unit of a rate and is linear in the 
scattering length a but depends in general on tempera¬ 
ture and other parameters (see below). 

Physically, Eq. describes independent oscillations 
of the two species in the absence of interactions. The 
eigenfrequencies of are given by the trap frequencies, 
TcJi and ±ci; 2 . Eq. 0 describes that each species intro¬ 
duces a Hartree potential for the other species. As we will 
discuss below, this contribution will shift the oscillation 
frequencies as long as the two species do not oscillate in 
parallel. W e ob tain withi n a local density approximation 
using Eqs. (|B9|) and (|B3|) from the appendix. 


5/23/2 .00 


7 = 


akBTuj{uj 2 m-Y m\ 
37rfi^ A' 2 ^red 


poo 

/ drr‘^gi{r)g2{r) 
Jo 


with 


fli(r) = Lii 


(8) 


( 9 ) 


where Lii is the poly logarithm of order ^ and yi is the 
chemical potential for particles of species i in the limit 
a ^ 0, see the discussion in appendix |B 1[ 

Damping, described by E(c(;), arises only to second or¬ 
der in the interaction strength. The total momentum is 
conserved during scattering processes, dtPi = —dtP 2 ^ 
which leads to the simple matrix structure 


with 


T = 


E(cj ^ 0 ) 


irh / 

M2kBT \2mred / 


/O 0 0 0 

0 0 0 0 
0 0 M 2 /M 1 -1 

Vo 0 -M 2 /Ml 1 


( 10 ) 


h n / 

n -1 O. 


= 1,2; ' 
a=l ,2 


d kiof 


X (5(Ae)(5(^)(Ak)g^ /ii/2i(l - /i2)(l - /22) (H) 
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to second order in the interaction strength using again 
the local density approximation, see Appendix |B 3[ Here 
fia are Fermi functions evaluated at the energy eia = 
and q = kn — ki 2 is the change 
of momentum of the first species, while Ak and Ae is the 
change of total momentum and energy, respectively. As 
the oscillation frequency is assumed to be much smaller 
than all Fermi energies, we have used the limit cj ^ 0. 
Furthermore, we ignore all frequency shifts to order o? 
(arising from the Kramers-Kronig partner of F). A more 
subtle issue is that our approach also neglects the cou¬ 
pling of the COM oscillations to other modes oscillating 
with frequency uji for a ^ 0. This is justified as, in the 
presence of interactions, these modes decay rapidly, but 
formally breaks down in the limit of vanishing interac¬ 
tions. As discussed in more detail in the supplement, 
this approximation gives rise to small, but nominally di¬ 
vergent extra contribution to 5](co’), which do, however, 
not affect our results. 

Finally, the equal-time correlation matrix Co in Eq. (§ 
is evaluated in Appendix |B 1 | To linear order in a, we 
obtain Co = Cq°^ + Cq^^ where 


Oq — 


n/{Miu:l) 0 0 0 

0 1/(M2 w|) 0 0 

0 0 Ml 0 

Vo 0 0 Ma/ 


( 12 ) 


and 


/ M2^2 


^( 1 ) _ 7 

° MiL0luJ2 


MiojI 

-1 

0 

0 


-1 

Miul 

M 2 UJ 2 

0 

0 


0 o\ 

0 0 
0 0 
0 0/ 


(13) 


where 7 (x a is given in Eqs. ii-®. 


II. ANALYTIC RESULTS 


Three different regimes have to be distinguished when 
discussing how the interactions affect the COM oscilla¬ 
tions, depending on which of the three quantities 6uj^ \j\ 
and F is largest. First, in the ballistic regime (Sou ^ 
| 7 |,r) interaction effects can approximately be ignored 
and the oscillations of the two species are almost inde¬ 
pendent. Second, in the frictionless drag regime (I 7 I ^ 
Jci;,r) one species drags the other by the interaction- 
induced Hartree potential. Finally, in the friction domi¬ 
nated drag regime (F ^ I 7 I) the two clouds are cou¬ 
pled by friction and only a hydrodynamic COM oscilla¬ 
tion with small effective damping survives. 

For a quantitative calculation we have evaluated the 
integrals in Eqs. @ and ( [TT| numerically, see section 
Jn] and Eig. In the limit of very low or very high 
temperature, also an analytic calculation is possible. Eor 



Figure 1. (Color online) Numerical results for the quantities 
7 (Eqs. (| 8 }-Q) and F (Eq. |TT] )) for different ratios of N2/N1. 
Dotted gray lines are analytic predictions for T <C cf,! and 
T ^ eF,i? see Eqs. The analytic formula for 7 in 

the limit T ^ 0 given in Eq. (14) is only exact for Ni = N 2 
and underestimates the value of 7 for N 2 /N 1 ^ 0 by a factor 
of 64 /(SStt) 0.58. All curves were calculated with mi = m 2 , 
^ = 0 . 1 , and are independent of the total particle number 
in the chosen units. 


7 : 


128 

357r^ 


~ 0.37 kF^idoJ 


87 r(fcBr )2 2 


(14) 


where kp^i is the Eermi momentum of species i in the 
center of the trap with Eermi energy ep^i = kpj{2mi) 
determined for T ^0. The analytic formulas have been 
computed in the limit ^ 0 and for mi = m 2 . While 
the prefactor of F is valid for arbitrary ratios of N 2 and 
Ni as long as N 2 < W, the prefactor for 7 is only exact 
for Ni = N 2 but increases by less than a factor of 2 when 
N 2 /N 1 is reduced, see Eig.[^ Surprisingly, the estimates 
given in Eq. (14) are even valid when the temperature 
is larger than the Eermi energy of the second species. If 
the temperature is larger than both Eermi energies, in 
contrast, the scattering rate F drops with 1/T while 7 
vanishes with 1 /T^/^, 


7 

F 


kF^iauj / kpauj / 

24v^ ~ 24v^ 

{kp^idY _ {kpaY ^F,i 
Ott nksT Ott h{kBTY 


(15) 


where kp = \/2nrik^T/h is the thermal wave vector. The 
prefactors for the high-temperature limit of both 7 and 
F are valid for arbitrary ratios of N 2 and Ni as long as 
N 2 < Ni. 

In both regimes, F can be identified with the single¬ 
particle scattering rate of a particle of species 2 in the 
center of the trap. In the high-temperature regime, this 
can be seen by rewriting F ^ cFVn^ni in terms of the 
sc attering cross section a ^ of ^ the typical velocity 'r’th ^ 
^ksT/m and the density of particles of species 1 in the 
center, ni ~ ~ ^^(m/T)^/^. 
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A. Ballistic regime 

For very sniall interactions the two species oscillate ap¬ 
proximately independently of each other. More precisely, 
we require that the strength 7 of the effective interac¬ 
tion potential and the single-particle scattering rate F 
are both smaller (in magnitude) than the difference of 
oscillation frequencies, \^\ 6uj and F <C While this 
regime is usually not realized experimentally at low tem¬ 
peratures (without tuning interactions close to zero), we 
discuss it here for completeness. Note that this regime is 
always reached in the limit of high temperatures as long 
as Suj 7 ^ 0 . 

The remaining weak interactions lead to a small shift 
of the respective oscillation frequencies relative to the 
trap frequencies and to a finite, but long, lifetime of the 
two oscillatory modes. The complex eigenfrequencies are 
given by the eigenvalues of — iS(cj), where the matrix 
of retarded susceptibilities, Eq. (§, has poles. We find 
for the eigenfrequencies in the ballistic regime. 



cjf/( 27 r) 


, ,ballistic 
‘^i 


M2u>2 7 . M2 r 

MiCOi 2 * Mi 2 


(16) 


where we evaluated both the frequency shift (real part) 
and the decay rate (imaginary part) to lowest order in 
the interaction strength a. Both the frequency shift and 
the decay rate are much smaller than Suj in the regime 
where Eq. (16) is valid. Eor low temperatures, T ^ 
one can use Eq. (14) to obtain for the frequency shift of 
the order of 


N 2 

^ kpiauj <C Suj 
Ni 


(17) 


while the decay rate of the oscillations is essentially given 
by the single-particle scattering rate. 


1 


N2 {ksTf 
Ni h^CF,! 


{kpci)^ Suj. 


(18) 


Eor high temperatures, T ^ ^f,i, the frequency shift 
drops faster than the decay rate and is therefore difficult 
to observe. 

In appendix I B 4| we show that in the ballistic regime the 
memory matrix formalism does not reproduce a dephas¬ 
ing of oscillations which gives rise to an extra effective 
decay rate linear in the scattering length a. This failure 
of the approach can be traced back to the fact that in the 
limit a ^ 0 an infinite set of further slow modes exists 
which we did not include into the set of slow modes , 


see appendix |B 4| for details. 

In a cold-atom experiment, one can directly observe 
the response of the clouds in real time. Erom the theory 
side, the real-time response can be obtained by Eourier 
transformation of the susceptibility, Eq. We consider 
the following setup: for time t < 0 a constant force is 
applied to the first species. Equivalently, we set in the 
Hamiltionian (Eq. ([^), = 0 and r5(t) = tq = roBx 

with ro > 0 for t < 0 (e^, is the unit vector in x direction). 


Figure 2. Response to a constant displacement rj(t < 0) = 
roBx of the trap potential for species 1 that is switched of sud¬ 
denly at time t = 0. The solid black (dashed blue) line shows 
the expectation value (Rf{t)) ((Rf (t))) of the center position 
of the first (second) atomic cloud, respectively, see Eq. (19). 
The calculations were done for Ni = A 2 , mi = ?ti 2 , ^ = 0.1 
and 7 and F as specified for each case. For high friction 
r (last plot), the oscillations of the two species synchronize 
quickly despite the finite difference Suj of the respective trap¬ 
ping frequencies, and the remaining COM oscillation decays 
only slowly on the time scale r/Suj‘^, see Eq. (25). 


The force is suddenly switched off, r^(t) = 0, for t > 0. 
In Eig. [^the expectation value 

.0 

{Ri{t)) = Miu^ro Xii{t - (19) 

J —OO 

is plotted as a function of time for both species, i = 1 , 2 . 
In Eig.[^ an example from the ballistic regime is shown. 
Due to the finite interactions the two modes couple and 
a beating pattern emerges which is characteristic for the 
superposition of the two frequencies uji and uj 2 . All os¬ 
cillations decay on a time scale set by l/F. 


B. Prictionless drag regime 

In the ballistic regime, the approximate symmetry 
which protects COM oscillations is of no relevance. This 
is different in cases where interactions are sufficiently 
strong so that the first species drags the second one either 
directly by the Hartree potential (frictionless drag) or by 
dissipative processes (friction dominated drag). In these 
regimes, the eigenmodes are characterized by a long-lived 
mode of COM oscillations, where both atomic clouds os¬ 
cillate in parallel, and a mode of relative oscillations, 
which decays more quickly. If all particles synchronize 
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their oscillation, then one can expect that the COM mode 
is approximately described by undamped oscillations of 
a rigid body of total mass Mtot = Mi + M 2 oscillating in 
an effective potential ^(MiCo’^ + M 2 Ci; 2 )r^. The oscillation 
frequency in this limit is given by 


UJ 


( 0 ) _ 
COM — 


IMiujI^M2UjI 


M 


tot 


( 20 ) 


We first consider the frictionless drag regime, which is 
reached when the frequency shift described by Eq. (16) 


becomes larger than the difference 5uj of the trapping 
frequencies, jyj ^ Jcj, and at the same time interaction 
effects between the two species are dominated by the ef¬ 
fective potential rather than scattering, i.e., jyl ^ T. For 
low T, this regime is obtained ioi ^ kpa ^ • 

In this frictionless drag regime we can use perturbation 
theory in 5uj to calculate the frequency shift and lifetime 
of the COM oscillations. We obtain 


<^COM ~ 


( 0 ) 

COM 


2 MIM 2 

V 7 7W 


ML 


( 21 ) 


As expected, wcoM large 7 ~ a, as the 

increasing drag effect causes the two atomic clouds to 
oscillate more and more in parallel despite the small dif¬ 
ference Suj of their trapping frequencies. Defining Auj by 
the shift relative to we obtain for low T 


Auj 


N 2 1 6uj‘^ 

Ni kpCL uj 


< Suj 


( 22 ) 


where we used again Eq. ( p^ . While the frequency shift 
is proportional to Xjkpa^ the lifetime turns out to be 
independent of the interaction strength in this regime. 


1 N2 ( Suj 

-^ — I — 

TCOM Ni \ UJ 


^ {ksTf 
h? eF,i 


Auj Suj 


(23) 


Note that both Auj and are proportional {Suj)‘^ as 

frequency shift and decay only arise from the small con¬ 
tributions violating the symmetry which approximately 
protecs COM oscillations. 

For completeness, we mention that the complex fre¬ 
quency of the mode of relative oscillations [27] is given 
by 


CJrel ~ ^COM 2Mi ^ 

This mode is damped by the single-particle relaxation 
time T and obtains a large frequency shift of the order of 
kpcioj for low T. As discussed above, the formula above 
ignores extra dephasing effects, see Appendix [C] 

In Fig. the real-time response is shown in the 
frictionless drag regime using again Eq. Due to 

the strong repulsive interactions the two clouds repel 
each other such that {Ri{t < 0 )) is larger than vq and 
{R^it < 0 )) is negative. After the external force has been 


switched off at t = 0 , the first cloud moves towards the 
center, first pushing the second cloud further away. Af¬ 
ter some time, the relative motion of the two clouds has 
decayed, the oscillations lock into each other and only 
the COM oscillations remain. The decay of the latter is 
given by the tiny rate r\j T{Suj/j)‘^, see Eq. due to 

the approximate symmetry. 


C. Friction dominated drag regime 


Experimentally, the most important regime is per¬ 
haps the hydrodynamic regime, where friction dominates, 
r ^ I7I, Suj. For ksT epi^ this condition is fulfilled for 

kpa > ^ \l ’ which is, e.g., re- 

alized with realistic experimental parameters of kp^ia ~ 
0 . 2 , Nik N 2 - 10 ®, 7 ~ 0 . 1 , and ksT « 0.3eF,i’ (the 
Fermi energies are given by eF,i = Note 

that this regime is always reached in the thermodynamic 
limit defined by ^ 00 , ^ 0 with ep^i = const. 

Furthermore, we demand as above that kpa I. 

The complex eigenfrequency of the COM mode is again 
obtained from perturbation theory in Suj and has the 
form 


<^COM 


,( 0 ) 

"^COM 


2 M 2 f Suj^^ .Suj^\ 


(25) 


Similar to the frictionless drag regime, interaction effects 
are suppressed for large T as the friction synchronizes 
the oscillations of the two atomic clouds. For low T we 
obtain the decay rate 


The frequency shift is in this regime much smaller than 
the decay rate, 

Acj < (27) 

^COM 


and therefore difficult to observe. For low T one obtains 

Fig. demonstrates how efficient a large friction is 
to lock the motion of the two clouds into each other on a 
time scale set by I/T. After this microscopic time-scale, 
only the center of mass oscillations remain, which decay 
very slowly on the time scale T /{Sujp , see Eq. (25). The 
motion of the two clouds is locked perfectly into each 
other. 


III. NUMERICAL RESULTS 

A. Protocols and possible experimental setups 

Depending on the setup of the cold-atom experiment, 
there exist various possibilities to access the different 
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A) increasing T B) increasing a C) increasing 6uj 



\Y\/6uj \Y\/6aj \y\/6uj 

Figure 3. (Color online) Depending on which one of the 
quantities Suj, I 7 I, and F is largest, the system is either in the 
ballistic regime (B), the frictionless drag regime (FLD) or the 
friction dominated drag regime (FDD). The arrows show the 
trajectories of the system in the parameter space when the 
temperature, the interaction strength a, or the difference Sou 
of the trapping frequencies are increased. All trajectories are 
calculated for Ni = N 2 = 10® and mi = m 2 . A) ^ =0.1 
(0.01), kF,ia = 0.13 (0.06), and ^7 = 0-05 ... 10 (0.07... 30) 
for the solid black (dashed purple) trajectory, respectively. B) 
^ = 0.1 ( 0 . 01 ), = 0.2 ( 0 . 1 ), and kF,iCi runs from 0 to 

0.025 (0.12) for the solid black (dashed purple) trajectory, 
respectively. C) kF,ia = 0.1 ( 0 . 02 ), = 0.2 ( 0 . 1 ), and 

^ = 0.03... 0.4 (0.002 ... 0.03) for the solid black (dashed 
purple) trajectory, respectively. 


physical regimes described in section [Ilj First, by chang¬ 
ing the cooling protocol, it is possible to access a broad 
range of temperatures. Second, by using an Feshbach 
resonance one can tune the scattering length. Third, if 
one has an experimental realization where the trapping 
potential of the two species can be varied independently, 
one can directly tune 6uj. In Fig. we show how each 
of these methods leads to a different trajectory in the 
parameter space spanned by j/Suj and V/duj. 

To illustrate the various regimes, we will plot in the 
following sections. Figs. and[^ the imaginary part 

of 


XcomM = (1,1,0,0) xM (1,1,0,0)^. (28) 

This describes the response of the center of mass to forces 
acting on both species simultaneously. Experimentally, 
the susceptibility as function of frequency can, e.g., be 
obtained by observing the real-time dynamics followed 
by a Fourier transformation. 


B. Increasing the temperature 


While 7 decreases monotonically as a function of tem¬ 
perature, F vanishes for both T ^ 0 and T ^ oc 
and has a maximum at ksT ^ see Fig. and 


Eqs. (14)-(15). Therefore, two scenarios are possible 


when one increases T while keeping all other parame¬ 
ters constant. If interactions are weak, kF^i\oi\ ^ ^ 
(solid black trajectory in Eig. [^), then the system is in 
the ballistic regime for low temperatures, may reach the 



ksT/ep,! 


lm[xcoM(^)] 

l{m 5u) ep 


80 


:60 


:40 


:20 

-0 


^ = 0.01; kp.ia = -0.06 ^ = 0.01; /cp.is = +0.06 

lm[xcoM(<^)] 

Itl50 
:100 
:50 
-^0 

0.1 1 10 0.1 1 10 
ksTlep,! ksTjer,! 




Figure 4. (Color online) Imaginary part of xcom(cc;), as 
defined in Eq. (28), as a function of temperature for attractive 
(left panels) and repulsive (right) interactions. The upper 
(lower) panels correspond to the solid (dashed) trajectories in 
Fig. 114, respectively. Dashed vertical lines indicate crossover 
temperatures where F = jyj or F = Suj, horizontal dotted lines 
are the analytical predictions of Eqs. ( p^ , ( [20| , and ( 21 ) using 
Eqs. (14)-|T^. For the top panels, the system evolves with 


increasing temperature from the ballistic (B) to the friction 
dominated drag (FDD) and back to the ballistic regime, while 
for the lower panels, the frictionless drag regime (FLD) is 


reached at low T. Parameters: Ni = A ^2 = 10®, mi = m 2 , ^ 
and kF,ia as stated above each plot. 



Figure 5. (Color online) Resonance frequencies cUres and de¬ 
cay rates l/r of the two eigenmodes of the system as a func¬ 
tion of temperature, calculated from the real and imaginary 
part of the eigenvalues of Q — zE, respectively, see Eq. © . In 
all four panels, the blue graph corresponds to the mode with 
longer life-time r. Dashed vertical lines indicate crossover 
temperatures between the ballistic (B), the frictionless drag 
(FLD), and the friction dominated drag (FDD) regime. The 
left (right) column corresponds to the top (bottom) panel in 
the left column of Fig. and to the solid (dashed) trajec¬ 
tory in Fig. [^, respectively. Parameters: iVi = A ^2 = 10®, 
mi = m2. 










































































1.0 


friction dominated drag regime at intermediate temper¬ 
atures ksT eF,i provided that kF,i\ci\ ^ 

returns to the ballistic regime for high temperatures. If, 
on the other hand kF^i\ci\ ^ ^ (dashed purple trajec¬ 
tory in Fig. IK). then the system is in the frictionless 
drag regime at low temperatures. Increasing the temper¬ 
ature to the order of the Fermi energy will typically drive 
the system into the friction dominated drag regime un¬ 
less kF,i\(i\ ^ hCd/eF,i ^ At high temperatures, 

the ballistic regime is always realized. 

For a quantitative analysis, we consider two concrete 
systems corresponding to the two trajectories in Fig. |3|\. 
In both cases, Ni = N 2 = 10^ and mi = m 2 . In the 
first system, ^ = 0.1 and kF,i\ci\ = 0.13, while in the 
second case we use ^ = 0.01 and kFiM = 0.06. We 
evaluate the integrals in Eqs. ^ and ( |11[ ) numerically 
for these two systems. Due to the spherical symmetry 
of the dispersion relation and the trapping potentials, 
the 12-dimensional integral in Eq. 0 can be reduced 
to a five-dimensional integral, which we evaluate using a 
Monte Carlo integration. 

The different regimes can clearly be identified in plots 
of Im[xcoM(<^)], Eq. ( [2^ , shown in Eig.[^ describing ex¬ 
citations of the COM motion. The vertical dashed lines 
in Eig.|^ correspond to the crossovers from one regime to 
the other, see Eig. |^, while the horizontal dotted lines 
give the analytical predictions for oscillation frequencies. 
The upper two (lower two) plots in Eig. [^correspond to 
the solid (dashed) line in Eig. On the left side, we 
consider attractive, on the right side repulsive interac¬ 
tions. 

The ballistic regime is characterized by the presence of 
two peaks: the two clouds oscillate independently with 
different frequencies. In contrast, a single peak located 
approximately at ujcom ~ ^ characterizes the two drag 
regimes where the oscillation of the two clouds synchro¬ 
nizes. A second (much broader) mode describing relative 
oscillations does not show up as for Xcom(<^) we only 
consider a situation where both clouds are displaced in 
the same direction (see Eig.j^for a plot of both resonance 
frequencies as a function of temperature). Note that for 
the chosen paramters, the system is not very deep in the 
ballistic regime for low T. This does not only lead to con¬ 
siderable shifts of the oscillation frequencies (see below) 
but also affects the weight of the two modes: the mode 
which is in frequency closer to clearly dominates. 

In the low-temperature ballistic regime, the interac¬ 
tions increase (decrease) the oscillation frequencies as 
the curvature of the potential increases (decreases) due 
to the attractive (repulsive) interaction with the other 
species, respectively. Interestingly, the effect is opposite 
for the drag-dominated regimes, best visible for the low- 
temperature regime in the lower two panels of Eig. [^ 
This higher-order effect, well described by our analytical 
formulas (( [^ and (25), drawn as dotted lines in Eig. [^, 
arises from level repulsion from the mode of relative os¬ 
cillations. 
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Figure 6. (Color online) Ini[xcoM(ct;)], Eq. ( [^ , as a func¬ 
tion of the scattering length a (corresponding to solid line in 
Fig. 1^). Dashed vertical lines indicate the scattering length 
where F = (5ct; and separate the ballistic (B) from the friction 
dominated drag regime (FDD). The dotted lines are analytic 
predictions of the eigenfrequencies based on Eqs. (16) and 
(25), where we used the low-temperature limit, Eq. (14), for 
the values of 7 and F. Parameters: Ai = A 2 = 10®, mi = m 2 , 
ksT = 0.2eF,i, ^ = 0.1. 


Fig. [^ shows the (real part of the) resonance frequen¬ 
cies and the decay rates 1 /r of both eigenmodes of the 
system for the two cases corresponding to the left panels 
in Fig.[^ The maxima of the decay rate of the long-lived 
mode (lower curve in the lower panels of Fig. trace the 
crossover from one regime to the next. The minimum in 
the friction dominated drag regime, where 1/tcom is pro¬ 
portional to the inverse of the single-particle scattering 
rate F, thereby arises from the maximum of F displayed 
in Fig.[^ While, in the ballistic regime, both modes have 
a long life time, in the drag regimes, only one long-lived 
mode remains and the decay rate of the mode of relative 
oscillations shoots up. 


C. Increasing the interaction strength 


As 7 oc a and F oc a^, the system evolves on a parabola 
in the parameter space of Fig. when the interaction 
strength is increased. While for weak interactions, the 
ballistic regime and for strong interactions the friction 
dominated drag regime is always realized, the frictionless 
drag regime is only reached if ^{ksT /< 0.03. 

Fig. ^ s hows numerical results for Im[xcoM (^]5 
c.f. Eq. p 8 [), for the solid black trajectory from Fig. |^. 
Dotted lines are again analytic results of the eigen¬ 
frequencies. The analytic prediction overestimates the 
slopes of the eigenfrequencies in the ballistic regime since 
it was made based on the T ^ 0 limit of 7 given in 


Eq. (14), while the actual value of 7 at ksT = {).2eF^i is 
by a factor of 0.62 smaller. 


D. Increasing the frequency difference 5uj 

Since 7 and F depend only weakly on Jcj for ^ <C 1, 
the trajectories for increasing 5uj in Fig. IP are almost 



















9 



lm[xcoM(t^)] 
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Figure 7. (Color online) Im[xcoM(co’)], Eq. (28), as 
a function of the frequency difference buj in the case 
(kBT!— 0.4 > 0.07 (solid black trajectory 
in Fig. [^). At the dashed vertical line, buj — F, separating 
the friction dominated drag regime (FDD) from the ballistic 
regime (B). The dotted lines are analytic predictions of the 
eigenfrequencies based on Eqs. ra and ( [ 2 ^ , where we used 
the low-temperature limit, Eq. (14), for the values of 7 and 
F. Parameters: N\ — — IO 7 mi = m 2 , kpa = ±0.1, 

ksT = 0.2eF. Here, cf and kp denote the Fermi energy and 
wave vector evaluated at (5ct; ^ 0 , respectively. 


straight lines crossing at the origin. For low temper¬ 
atures and small Su, the friction dominated (friction¬ 
less) drag regime is realized if {kBT/€F,i)‘^ kF,i\ci\ 
is larger (smaller) than 0.07, respectively. For large Su 
(and weak interactions), the system enters the ballistic 
regime. Fig. shows numerical results for Im[xcoM (<^)]5 
c.f. Eq. ([^), corresponding to the solid black trajec¬ 
tory in Fig. Dotted lines are analytic predictions of 
the eigenfrequencies based on the low-temperature limit, 
Eq. (ig. 


IV. CONCLUSIONS 


but similar mass. Alternatively, one can also investi¬ 
gate, e.g., two spin species with the same mass but 
slightly different harmonic confinement. The mass dif¬ 
ference and/or difference in the strength of the parabolic 
potential breaks a dynamical symmetry which otherwise 
protects the center-of-mass oscillations from decay. 

The interactions of the two species synchronizes the 
motion of the two clouds and thereby leads to a partial 
restoration of the dynamical symmetry: the interacting 
liquid can approximately be viewed as having a single 
average mass and oscillating in a single average potential. 
As a consequence, the decay rate of the center-of-mass 
oscillations is strongly reduced and of the order of , 
where Su is the difference of the trapping frequencies and 
F the scattering-rate of the two species. Compared to 
other hydrodynamic modes (which can also have decay 
rates proprotional to the inverse of F) one obtains an 
extra reduction by the factor {Su/u^. 

As all other modes have much faster decay rates, the 
approximate symmetry leads to an almost perfect drag of 
the two clouds: the center-of-masses for each of the two 
species follow each other after a few scattering times. 

For future investigations two directions are especially 
interesting: First, one can study the highly non-linear 
regime, where, for example, initially one species is sep¬ 
arated far from the second one and one can study the 
evolution of the center-of-mass oscillations after the two 
clouds have violently crashed into each other in a setup 
similar to the one studied by the Zwierlein group [38]. 
Second, one can investigate the interplay of superfluidity 
and the approximate symmetry, which is of direct rel¬ 
evance for the experiments of the Salomon group m- 
Here, in the center of the cloud and for small relative ve¬ 
locities, the superfiuid components move without friction 
and only the normal components can scatter from each 
other. 


The presence of approximate symmetries leads to a 
slow equilibration of a perturbed system. We suggest 
that this physics can be studied with high precision ex¬ 
perimentally by investigating the center of mass oscil¬ 
lations of two species of ultracold atoms with different 
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Appendix A: Memory matrix formalism 

For a general set of observables time-dependent 
forces fn{t) on each observable are described by a contri¬ 
bution iLext = to the Hamiltonian. The 

response of some observable Am to the external forces 
fn{t) is described, to linear order in /^, by the matrix 
of retarded susceptibilities Xmn{^)i as defined in Eq. 
via the relation 

{Am{uj)) = 2 tt5{uj){A^) 

eq. 4“ E 

n 

where the Fourier transform of the external forces (and 
accordingly of the observable Am) is defined by fn{^) = 


fit)dt and (•) denotes the expectation value in 
the perturbed system while (•)eq. is the equilibrium ex¬ 
pectation value for /n(t) = 0 Vt. 

To be specific, as described in the main text we study 
the dynamics of the (x-components of the) COM coordi¬ 
nates of two atomic clouds. Thus, An = {Rf, Pf ^ Pf) 
as defined in Eq. The Hamiltonian, Eq. § , contains 
forces on the COM position coordinates Rf, i = 1,2 of 
the two atomic clouds, which are given by m^co’|r^’^(t). 
We do not consider forces on the COM momenta Pf, 
but nevertheless include Pf in the set of operators An 
since (i) the momenta can be observed in time-of-fiight 
measurements and (ii) we expect excitations of the COM 
momenta to be long-lived in the regime ^ ^ 1 and the 
applied memory-matrix formalism requires a separation 
of time scales where the operators An span the subspace 
of all slowly relaxing local observables. 

We calculate the matrix of retarded susceptibilities 
Xmni^)i E g. Q , by means of the memory-matrix for¬ 
malism [29143 1] 7 In the following, we briefly review the 
central results of this technique. 

A scalar product in the space of quantum-mechanical 
operators is defined by 

r/3 

{A\B):= dX{A^B{ihX)U^,-(3{AXci.{B)eci. (A2) 
Jo 

where (3 = l/iksT) is the inverse temperature and 
B{ihX) = is the operator in the Heisenberg 

picture. Instead of calculating the matrix of retarded 
susceptibilities Xmni^) directly, it is easier to first de¬ 
rive an expression for the matrix of retarded correlation 
functions Cmni^) defined by 

nOO 

Cmn{oj)= (A3) 

Jo 

It is easy to show that Xmni^) and Cmni^) are related 
via 

Xmn(^) — tujCmni^) 4~ (^o)mn (-'^d) 

where the entries of the equal-time correlation matrix Cq 
are defined by 

(t^o) mn — {Am\An). (^5) 

Time evolution of an operator is described by A{t) = 
e^^^A with the Liouville (super-)operator L = ^[P, •]• 
Thus, Cmni^^) is given by 

pOO 

Cmnico)= 

Jo 

= - A)“Mn) (A6) 

for Im(cc;) > 0. 

The operators An span a subspace of the space of 
quantum-mechanical operators. We define the projection 


11 


(super-)operator V (Q) onto (away from) this subspace 
by 


V=l-Q = Y,\Am){C^^)mn{An\. (A7) 

m,n 

Inserting L = LQ-{-LV into Eq. (|A 6 | ) and following some 
simple algebraic manipulations |29j one arrives at a ma¬ 
trix equation for the retarded correlation functions, 

C(w) =i(w-fi + iI](w))“^Co (A 8 ) 

where 

^mn = i'^{Am\As){CQ'^)sn (A9) 

S 

Smn(w) = iJ2i^rn\Q{0J - LQ)-^\As){Co^)sn (AlO) 


The matrix O describes the evolution of the observables 
An{t) if there was no coupling to any other degrees of 
freedom (i.e., if L would commute with all An). Effects 
due to the coupling of the An modes to other modes are 
encoded in the memory matrix E(ci;). 


Appendix B: Evaluation of the matrices Co, 11, and 

S(c^) 


In this s ection we eval uate t he matrices Co, H and 
I](w), Eqs. ( |A5| ),( A9), and (AlO), for the model described 
by Eq. (§. All calculations are done perturbatively for 
small interaction strength a. Scalar products are calcu¬ 
lated in the local density approximation, which is valid 
for Ni,N 2 > 1 . 


1. Equal-time correlation matrix Co 


T he e qual-time correlation matrix Co is defined in 
Eq. (A5). Due to their different signature under time 
reversal, the position and momentum operators have 
vanishing overlap, {Rf\P^) = 0. We expand Co ~ 

Cq^^ -h + O(a^) for small a. Without interactions, 
a = 0 , the two species decouple from each other, result¬ 
ing in a diagonal matrix structure of We obtain, 

within a local density approximation. 


(Bl) 

0 0 0 M 2 / 

where Mi = Nirrii. 

In local density approximation, the momentum- 
momentum components of Co are not changed by inter¬ 
actions to first order in a. Interactions only affect the 


Oq — 


/l/{Miujf) 0 0 0 

0 1 /(M 2 Cj|) 0 0 

0 0 Ml 0 


position-position components and we obtain 


f- 






M2UjI 

1 

0 

0 


M 2 UJI 

0 

0 


0 o\ 

0 0 
0 0 
0 0/ 


(B2) 


where 




1 akBT{jTi\m2)^l^ 


f 


A’iA"2mred 
gi{r) = Lii . 


drr‘^gi{r)g2{r) 


(B3) 


Here, Lii is the poly logarithm of order | and is the 
chemical potential for particles of species i in the limit 
a ^ 0. To arrive at the diagonal matrix elements of 
given in Eq. (B2) one has to take into account that 


the actual chemical potentials depend on the interaction 
strength a. Eq. (B2) gives the result for fixed particle 


numbers, which was derived using the relation 



d{Rm^) 

dNi 

da 

Ni 

dNi da 


Mi 

(B4) 

where the notation denotes that x is kept constant in 
the derivative. Eor the off-diagonal elements of Co, cor¬ 
rections due to the dependency of the chemical potentials 
on the interaction strength are of higher order in a. 


2. Eigenfrequency matrix Q 


We expand D + O(a^) for small a. The 

temporal derivatives A m = ^ [M, Am] that appear on the 
right-hand side of Eq. ( |A9[ ) are given by 


/Mi 

Pf = -MiUjjRf 


(B5) 


^(12) px 
^int 


=: P- 


,trap 


JDX 

M,int 


Setting = 0 in Eq. (B5) and inserting into Eq. (A9) 
using Eq. (A5) one arrives directly at the leading-order 
contribution to the eigenfrequency matrix. 


/ 0 0 i/Mi 0 \ 

0 0 0 i/M2 I 

-iMiuol 0 0 0 

V 0 —iM2UJ2 0 0 / 


(B 6 ) 


Note that, once we set Pf-mt = 0 in Eq. ( |B5[ ), all scalar 


products that appear in the evaluation of are can¬ 
celed exactly by the factor on the right-hand side 
of Eq. (A9). Thus, all corrections to due to inter¬ 
actions originate from the term i^ Eq. (B5). Since 


due to different signature under time re¬ 
versal, interactions only change the lower left 2 x 2 corner 
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(X qiq2 X 



0'?r<C>-r'0 


Figu re 8. The two types of diagrams contributing to 
Eq. ( |B12| ). Solid (dashed) lines represent fermions (inter¬ 
species interactions), respectively. The second diagram would 
vanish in a homogeneous system due to momentum conserva¬ 
tion at the vertices. In a harmonic trap, however, E^“^^ does 
not vanish and has poles at the trap frequencies (for vanishing 
intra-species interactions). 


of the matrix Q. We get, to leading order in a, 

(A",intl^ 2 ) = -(A^intl^l) = -MlLol (B7) 

Where is given in Eq. ( |B3| ). Inserting Eqs. ( |B7| ) 

into Eq. ( |A9[ ) yields the first order correction to the eigen- 
frequency matrix 


=i{Rf\R2) 


-M 2 U 0 I 
\ M 2 UO 2 


0 0 0\ 

0 0 0 I 

Miul 0 0 po 

-MiujI 0 0 / 


(B 8 ) 


Since the factor {RflR^) is already first order in a we 
may approximate the matrix Cq by (Eq. ig). This 
leads to 


a we e valuat e all scalar products on the right-hand side 
of Eq. (AlO) with respect to the non-interacting system 
and describe time evolution by the non-interacting Li- 
ouvillian Lq = [Hq^ •]. Using further the fact that Lq 
commutes with Q and that Cq is diagonal (to lowest or¬ 
der in a) one finds that only the P, P-components of E(co’) 
have non-vanishing values given by 


Epcc px(uj) = I -^ ^ ^ , - - - 


Inserting Q = 1 —V leads to 


px px 


(cj) — px (cj) + E 


(b) 

px 


,(w) 


(BIO) 


(Bll) 


with 


-,(a) 

px p: 


(w) = 


* iP!\P!) 

H + (w) 


px p. 


^{uj) = -i 


.iPi 


i,int 


\r{x - 


(PfIPf) 


(B12) 

(B13) 


Two types of diagrams, shown in Fig. contribute to 
E;(^), which we denote by E^^^^ and E^^^ in the follow¬ 
ing. E^^^^ describes how the scattering of quasi parti¬ 
cles leads to momentum transfer from one species to the 
other. Evaluating the diagram within the local density 
approximation (i.e., by approximating the system locally 
by a homogeneous one) results in Eqs. (lO)-(ll) of the 
main text. 


4. Singular contributions to the memory matrix 




/ 0 0 0 0 \ 
0 0 0 0 
1-10 0 
V-1 1 0 0/ 


(B9) 


Eqs. (B 6 ), (B9), and (B3) describe our result for the lead¬ 
ing and next-to-leading order contribution to the eigen- 
frequency matrix Q. 


3. Memory matrix E(c< 


T he def inition of the memory matrix T^{uj) is given in 
Eq. (AlO). Note that one may insert an additional pro¬ 


jection operator Q to the le ft of the vector jA^) in the 
right-hand s ide of Eq. (AlO) without changing its value. 
Using Eqs. (|B5[) and the fact that Q projects onto the 


subspace of obser vable s orthogonal to the space spanned 
by \An) (see Eq. (|^) we find Q|Pf) = 0 = Qli^trap)- 
Therefore, all contributions to E(co’) come from terms 
quadratic in Pfi^t thus at least of second order in the 
interaction strength a. Neglecting higher order terms in 


The discussion of the second diagram in Fig. § 
and of requires some more care. While both contri¬ 
butions vanish within the local density approximation, 
they have both a divergent contribution for uj = uji \i 
evaluated exactly (as long as intra-species interactions 
are absent, see below). We will argue that the diver¬ 
gencies cancel exaetly if the full space of slow operator 
is considered. This is the reason why we ignore these 
contributions for our analysis despite the fact that the 
cancellation is only partial for the four modes considered 
by us. 

First, we note that E^^^^ has a simple physical interpre¬ 
tation. It describes that the momentum decays because, 
e.g., species 1 is affected by a single-particle Hartree po¬ 
tential VA(r) = |^^( 4 ^ 2 (^)^ 2 (i*)), giving rise to an ex¬ 
tra force Fhi = — f 4^|(r)4^i(r) d^r contributing 

to dfPi. We obtain 


y,(a 2 ) 
2^px px 

% ■> J 


(w) + Spl px{to) 

i ’ j 


- Lo)-^\FHj) 

{PjpPj^ 

(B14) 
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In cases where we can approximate Vh by a parabola, 
¥h is proportional to Ri and does not contribute to H as 
QFfji = 0 . Similar, if we include in our list of slow modes 
An all operators which oscillate with frequency then 

by construction Y^^px^px (uj) + Ypl px (uj) is non-singular 
for uj Fuji. The fact that for our choice of slow modes 
an extra divergency remains, arises from a peculiar prop¬ 
erty of the harmonic oscillator. As all single-particle en¬ 
ergy levels are equally spaced, there is an infinite number 
of hermitian operators oscillating with frequencies Fuji in 
the non-interacting limit, 

= j ■^\{v)r,hT-^i{v)Sr (B 15 ) 

= -ij (B 16 ) 


where hi = is the single particle 

Hamiltonian of species i. 

In our analysis we have (i) only included the operators 
with m = 0 and (ii) neglected the divergent contributions 
discussed above. It is therefore important to ask to what 
extent our results are modified when further terms with 
m > 0 are included. First, the accuracy of the result will 
increase as by construction the neglected terms become 
smaller and smaller. It is important to note, that the 
more complicated operators Am^ij with m > 0 are not 
protected by any approximate symmetry. Therefore their 
decay rate is not suppressed by factors of { 5 uj/ujY and 
they will neither qualitatively nor quantitatively influ¬ 
ence the final results in the limit where interactions be¬ 
come important (the hydrodynamic friction dominated 
drag regime). Furthermore, in situations where intra¬ 
species interactions are present, these lead to a decay of 
with m > 0 but do not affect Am^ij with m = 0 . 
Therefore, if the decay rate due to intra-species interac¬ 
tions is sufficiently high, our results are again fully valid 
in all regimes considered. Techni cally, this is reflected by 
the fact that the omitted terms (BI 4 ) are non-divergent 
if intraspecies interactions are included in Lq. 

In the ballistic limit where all friction can be ignored, 
however, our results presented in the main text miss a 
physically important effect: the decay of oscillations by 
dephasing (rather than decay by friction considered by 
us). As we will show in the following section, this leads 
in the ballistic regime to a decay rate which is linear in 
the scattering length a and is not covered in memory 
matrix approximations which neglect the modes Am^ij 
with m > 0 and predict decay rates proportional to a^. 


Appendix C: Decay by dephasing: a toy model 

In this section, we discuss the decay of COM oscil¬ 
lations for a simple toy model where non-interacting 
fermions scatter from a weak, smooth, and time- 
independent potential V{r). We will use the calcula¬ 
tion to show that (i) the memory matrix approach cor¬ 


rectly describes the average shift of frequencies due to 
the Hartree potentials but (ii) fails to reproduce the cor¬ 
rect lifetime in the ballistic regime due to the problems 
discussed in section [B 4 ] above. 

In the simplified model considered here, the atoms of 
the first species do not move while those of the sec¬ 
ond species oscillate in a potential given by the har¬ 
monic trap ^m2002r‘^ plus the static Hartree potential 

V{r) = ^ (ni(r)), where (ni(r)) is the expecta¬ 

tion value of the density of species I in equilibrium. 
Formally, we consider the limit rui oo, miujf = 
1712002 , Ni = ^2(7711/7712)^/^, kF,2a (7772/7771)^/^, such 
that the two clouds have the same shape (fii(r)) 
(mi/m2)^/^(n2(r)). 

The imaginary part of the retarded susceptibility for 
the COM position of the second species is given by the 
Kubo formula, 

^{fOa) - /(£«')) Ka'|4|a)px 

^ a,a' 

X 6 {tuJ - (Ca' - Ca)), (Cl) 


where / is the fermi function and ja) are single-par ticle 
eigenstates with energies Cq,. We evaluate Eq. (Cl) 
perturbatively for small V. As the energy levels 


= hjoj2{n + |) of the unperturbed three-dimensional 
isotropic harmonic oscillator are degenerate, one has 
to diagonalize the matrix {a'\V\a) for each 77-subspace. 
Since V{r) oc ( 77 i(r)) is spherically symmetric, this is 
done by the states ja) = 1 77/777), where / and m are the 
quantum numbers of angular momentum. To linear order 
in V, the eigenenergies of these states are independent of 
777 and given by Cn^i = + ( 77 /|E| 77 /), which we evalu¬ 

ate numerically by a one-dimensional integration in the 
radial direction. As we are only interested in the be¬ 
havior of lm[xR^,R^{oo)] for uo close to UO2 and to linear 
order in V, it is sufficient to calculate the matrix ele¬ 
ments {n'l'm'\rx\'alm) to order which leads to the 
selection rules 77' = 77±I, /' = /±I, and m' = m (with 
quantization axis in the x direction). We obtain 


+ 1,1 ±l,m\rx\n,l,m)f = 

m 

_Vh{2l + 1){21 + 1 ± 1){21 + 1 ± 2 )( 2 n + 5±2l±l) 

~ v^miu^24(2Z ± 1){21 + 2 ± 1) ' 

(C 2 ) 

In Fig. 1^ we show the resulting Ini[xijj,_Rj (w)] for x 
close to 002 ^ The (5-peak one would obtain for E = 0 at 
cj = Cc ;2 is shifted linearly in the scattering length a and 
also broadened linearly in a. Note that the peak shape 
is not Lorentzian. The broading is not caused by inelas¬ 
tic scattering but arises instead from a simple dephasing 
effect: the frequency shifts linear in V affect the energies 
of different eigenstates in a different way. The dephasing 
linear in V is not covered by the version of the memory 
matrix approach used by us, which does not take into 
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Figure 9. Imaginary part of the retarded susceptibility for a 
non-interacting gas in a harmonic trap disturbed by a small 
extra potential V(r) — ^(hi(r)), see Eq. (Cl), where the 

(5-function is approximate<f by a Gaussian with tiny standard 
deviation — 0.01|Act;|. The das hed line is the average shift 
Ao; of the peak position, Eq. ( |C3|) , correctly predicted by the 
memory matrix method, Eq. ( |C5| ). Parameters: /csT/eF,2 = 
0 .1, N 2 — 10^ in the limit mil m 2 00, miuJl — m 2 UJ 2 ^ 

Ni = A2(7711/7712)^^^, kF , 2 a <C (7712/7711)^/^. 


with the normalization C = ^\m[xR^,R^{w)] = 

1/(4A^2^2<^2) and 5 < UJ 2 chosen such that only the 
weig ht o f the peak close to co ’2 is captured. Inserting 
Eq. ( |C1[ ), we find to linear order in 

^ - /(4°4) I(a'l4|a)l^x 

ol^ol' 

X {{a\V\a') - {a\V\a)) (C4) 

where d|. is the ladder operator of the harmonic oscillator 
in X direction. 

Applying our version of the memory matrix, we obtain 
for the frequency shift to linear order in V 

a,a' 

(C5) 


account higher modes Am,i,j with 777 > 0, and ignores di¬ 
vergent terms in T>{uj) arising in a treatment beyond the 
local density approximation. This is the main result of 
this section. Note that this dephasing only affects the de¬ 
cay rates in the ballistic regime, Eq. (16) and the decay 


rate of the relative motion of the two species in the drag 
regimes, Eq. (24). It is irrelevant in the friction domi¬ 
nated drag regime, where the inelastic scattering rate E 
is much larger than the dephasing rate. Furthermore, the 
dephasing effects are expected to be strongly reduced by 
the factor for the COM mode in the frictionless 


drag regime, cf. Eq. (21), as for the synchronized motion 


of the two species the Hartree potential cancels to leading 
order. 

In the following, we will show that our memory ma¬ 
trix approach does, however, correctly predict the aver¬ 
age frequency shift, see dashed line in Fig.[^ The average 
frequency shift is defined by 


1 pUJ2-\-S 7 

Ao; = 77 / X- {uj - uj2)lm[xR-,R-{uj)] (C3) 
C Juj2-6 27r 


where Ff = —f d^r is the force aris¬ 

ing from the Hartree potential, and we used ^ = 
y^2m2002lh[dx^ V] in the last e qual ity. We factorize 
the last matrix element in Eq. ( |C5[ ) by inserting 1 = 
|<a)((a| between the operators dx and and write it 
explicitly in the eigenstates \nlm) of the harmonic oscil¬ 
lator, 

{a\[dx^V]\a') = (^{nlrn\dx\nlm){nlm\V\nd'rn') 

h,l,fri 

— {nlm\V\hlm) {hlm\dx\n'dm')). 

(C6) 


Inserting Eq. ( |C6| ) into Eq. (C5), using the properties 
{nlm\V\nlm) oc (5^[^rn,m for the spherically symmetric 
potential V and {nlm\dx\film) oc ^n+i,n? reproduces ex¬ 


actly the average frequency shift, Eq. (|C4|), derived from 
the direct calculation of lm[xR^,R^{uj) 


















